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Collective excitations both for spin- and charge-channels 
are investigated in incommensurate spin density wave (or 
stripe) states on two-dimensional Hubbard model. By random 
phase approximation, the dynamical susceptibility w ) is 
calculated for full range of (q, a>) with including all higher 
harmonics components. An intricate landscape of the spec- 
tra in x(Q) w ) is obtained. We discuss the anisotropy of the 
dispersion cones for spin wave excitations, and for the phason 
excitation related to the motion of the stripe line. Inelastic 
neutron experiments on Cr and its alloys and stripe states of 
underdoped cuprates are proposed. 

PACS numbers: 75.30.Fv, 75.10.Lp, 72.15.Nj 



Recently a remarkable scries of elastic neutron exper- 
iments has been performed on I^-xSr^CuC^ 0. It re- 
veals static incommensurate spin density wave (ISDW) 
structure in underdoped regions. The superconducting 
samples show that the modulation vector Q are char- 
acterized by Q = 27r(i, A ± S) or 27r(A_ ± 5, A) in mo- 
mentum space. The incommensurate modulation runs 
vertically in the a- or &-axis of the CuC>2 plane. The 
static stripe structures have been observed also in su- 
perconducting (La, Nd)2-xSr x Cu04 p| and insulating 
La2- :E Sr ;c Ni04 M. Since the establishment of static or- 
ders on these systems, the study of spin and charge dy- 
namics; x(q, a;) has been just started experimentally. 

Prior to these studies on cuprates, ISDW is well known 
in itinerant electron systems, such as a typical example 
Cr and its alloys for long time. While its static proper- 
ties are fairly well understood, its dynamical properties 
remain largely unexplored QJj| . A part of reasons stems 
from the fact that the theoretical description is under- 
developed, thus the interplay between theory and exper- 
iment was not satisfactory ||. The intensive efforts for 
extracting spin and charge excitations in Cr and its al- 
loys and also high T c cuprates are strongly motivated by 
a hope that knowing the collective excitations or fluctu- 
ation spectrum in strongly correlated systems may lead 
us to a clue to understanding the mechanism of high T c 
superconductivity. 

There exist a lot of theoretical works on dynamical 
properties, beginning from a seminal paper by Fedders 
and Martin to recent work by Fishman and Liu |l. 
The former first derives spin wave mode in the transverse 
spin excitation of itinerant electron systems. The latter 
investigates various transverse and longitudinal modes at 
the Q position, based on a one-dimensional (ID) model 



within the random phase approximation (RPA). Their 
theory takes account of only fundamental order parame- 
ter Aq , neglecting higher harmonics A2Q , A3Q , • • • as- 
sociated with the incommensurability S. By extending 
their work |jj , we calculate the dynamical spin and charge 
susceptibilities xCli w ) m the full range of the two dimen- 
sional (2D) wave number q for entire Brillouin zone and 
the excitation energy u> up to the band width. We take 
account of all the possible higher harmonics. This will 
turn out to be extremely crucial in correctly evaluating 
these quantities. The multi-dimensional calculation here 
allows us to extract a wealth of information on stripe mo- 
tions such as translation, or meandering, etc. and on the 
anisotropy of excitation cones. This kind of calculation 
has not been done before to our knowledge. 

We start with the Hubbard model on 2D square lat- 
tice: H = ' V,., „ <",.,/',.„ + UY.i nifrni. To calcu- 
late %(q, u>), we first set up the incommensurate SDW 
ground state. Assuming a periodic spin and associ- 
ated charge orderings, we introduce the order parame- 



ter (n it a) = J2i e 



iiQ-i-i 



HQ,tr) with 5 = 1/N. In AT-site 



periodic case, the Brillouin zone is reduced to l/Af-area. 
and energy dispersion is split to N bands. We write 
k = k + mQ (m = 0, 1, • • ■ , N — 1), where k is 
restricted within the reduced Brillouin zone. Then the 
Hamiltonian is reduced to 



H ~ Ck +mQ,o'(^ko,cr)m,m'Ck +m'Q,<T 
k ,<T m,n 



^ko,cr,a7k ,£r,ct7k ,a-,a- 



(1) 



ko,c,a 

The N x N Hamiltonian matrix 

(#k , o)m,m> = e(k + mQ,)S m>m : + m')Q, — cr) IS 

diagonalized by a unitary transformation Ck +mQ,tr = 
X] a Wk ,cr,a,m7k,er,a- The calculation is iterated until all 
the order parameters satisfy the self-consistent condition 

Here,AT fe = E k0>ro l. 

We construct the thermal Green function as 

/ / • v V- "k , ff , a (rK ^ a (r') 
g a (r, r , iu n ) = ^ — ^r 2 ^ , ( 2 ) 



k ,Q 
-1/2- 



iu! n — i?k 



a i(k +mQ)-r i 



and 



With Uko.a.aOi) = N k 1 E m e '^ UT ""^' lw k ,<T, Q ,m, 

evaluate the dynamical susceptibilities; the spin longitu- 
dinal mode x zz (q, lo) = ((S z ; <5 , z ))q,u J , the transverse one 
Xxx{<\,u) = ({SxiSx))^, and the charge susceptibility 
Xnn(q, — ({n\ n)}q.cj. The Fourier transformation of 
Eq. (0) is given by 



1 



g a (k + ZiQ, k + Z 2 Q, iw n ) 



E 



ri,r 2 

iw„ — £/k.cr,a 



(3) 



where /i, Z 2 = 0, 1, ■ ■ ■ , N — 1. In the presence of the 
order parameter {niQ l0 -}, the incoming- and outgoing- 
momentum of g a can differ by ZQ = (Zi — Z 2 )Q. 
Then g a {k + ZiQ,k + Z 2 Q,kj rl ) is an N x N matrix 
with indexes Zi and 1%. The bare susceptibility is 
given by Xo CT '( r i, r 2,i^«) = -^E^^fn-^.iwn + 
ifi„)<7 CT '(r2,ri,kj n ). Its Fourier transformation Xo CT (k + 
ZiQ, k + Z 2 Q, is also N x N matrix 0]. We use the 
analytic continuation kj n — > cj + i?7. Typically, we use 
7? = 0.001£ in our numerical calculation. 

The RPA equation for XSfSi = ((^ti is written as 



+U S T / xi L ( r i > r 2 . w)xs T s ; (r 2 , r 3 , u>) . 
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(4) 



After Fourier transformation to k-space, Eq. (0) is re- 
duced to a matrix equation of Xo and XStSi • By solving 
it, we obtain Xxx(k + ZiQ, k + Z2Q, to). In a similar man- 
ner, we calculate {{n-\;n-\)) and ((n-\;ni)), and obtain 
Xzz( k + 'xQ>k + Z 2 Q,w) and x«n(k + Z x Q,k + Z 2 Q,w). 
The neutron scattering experiments observe the imag- 
inary part of the dynamical susceptibility x"(q> w ) — 
Imx(q, q, u>). Since the signal is observed as spatial av- 
erage, we consider the diagonal part of x(k + ZiQ, k + 
Z2Q, w). 

According to standard linear response theory, the 
spatio-temporal oscillation pattern of a collective mode 
can be analyzed by x- I n the presence of an infinites- 
imal external field Zio'(k',w) coupled to an operator 
0'{= S x , S z , n), the response of the operator O is 
given by S(O(k + hQ,w)) = Xoo'iK + ZiQ,k + 
Z 2 Q,w)/io'(ko + Z 2 Q,w). When the external field is a 
plane wave; ho<{v,t) — foo ,e ^ q ' r_ with a small am- 
plitude ho> , the response is given by 



5{0(r, t)) = £ Xoo> (q + IQ, q, co)e^ +l ^h 
1 



h Q ,e 



— iujt 



(5) 



We consider the vertical stripe case for U = At and 
the hole density — 1/8 as a representative case (The 
enrgy is scaled by t from now on). As we do not in- 
clude the nearest neighbor hopping t' , the lowest energy 
ground state is an SDW-gapped insulator with 6 — 
i.e. N — 16. The detailed ground state properties 
are reported previously pl]Jl2|. For our parameters, 



the single particle SDW gap E g = 0.41. In the ISDW 
state, the spatial profile is characterized by a distorted 



sinusoidal, or soliton form with a midgap band. The 
higher harmonics are determined as \Miq/Mq\ =0.08 
(Z = 3), 0.005 (Z = 5), 2.3 x I0~ 3 (Z = 7), where 
Miq = (jfyQ.-f) — (jizq,j.)- In the limit of the half- filling 
(N — > 00), the ratio |M( 2n+1 )Q/Mq| increases and ap- 
proaches (2n + since the profile of the spin struc- 
ture approaches square wave form jTl],[l2| . The structure 
factor has spots at (2n + 1)Q in the spin structure, and 
at 2nQ in the charge structure. These spots are observed 
by the elastic neutron scattering in Cr and its alloys 0. 
Their positions in the momentum space are shown in Fig. 

0- 

Let us start with the excitation of the spin transverse 
mode. In Fig. ||, we show x"x(q> w ) along paths A and 
B of Fig. El The gapless spin wave modes emanate 
not only Q, but also from 3Q and other odd harmon- 
ics. The ridge of x"(q> w ) shows singularity reflecting 
the dispersion of the collective mode. These modes at 
(2n + 1)Q have an identical dispersion relation, because 
every (2n + l)Q-modes couple each other in the RPA 
equation (|j). The same dispersion pattern appears in 
each reduced Brillouin zone. But their intensities are 
different. With increasing Z, the intensity is decreased 
as follows, x"x('Q:0-025i)/x"x(Q,0.025t) = 6.5 x I0~ 3 
(Z = 3), 2.7 x 10~ 5 (Z = 5), 5.8 x 10~ 6 (Z = 7). It is found 
that the intensity ratio is obeyed 

2 



X^((2n + l)Q,^->0) 



M, 



(2n+l)Q 



Mi 



Q 



(6) 



With approaching the half-filling (N —> 00) the above 
ratio is expected to increase and approach (2n + 1)~ 2 . 
Then, the higher harmonics spot at (2n+ 1)Q may have 
enough intensity to be observed near half-fillings such as 
in Cr or underdoped cuprates. 

We analyze the oscillation pattern by Eq. (0). The 
response of 5{S x (r, t)) shows the same spin pattern as S z 
of the ground state at to ~ for (2n+ 1)Q. It means that 
the ground state spin structure is rotated as it is without 
modulation, since it is a Goldstone mode. The external 
field of the wave number (2n + 1)Q couples to M(^ n +i)Q 
and makes the ground state spin structure rotate. Then, 
we can conclude that the spin transverse mode is a spin 
wave. With increasing ui along the dispersion curve, the 
spin wave oscillation shows the deviation from the spin 
pattern of the ground state, reflecting the wave number 
of the external field. 

The reconnection of the dispersion curve occurs at the 
reduced Brillouin zone boundary. Then, instead of the 
simple intersection of two dispersions, the small gap ap- 
pears at ui ~ 0.3i in Fig. |^. With increasing cj, the 
intensity of (<!?<*>) decreases as 1/uo along the disper- 
sion curve, except for the weakened intensity at the gap 
position. For ui > E gi there exists other modes reflecting 
the fluctuation of the magnetic moment amplitude. It is 
a character of itinerant magnets and absent in localized 
spin magnets. 
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Figure [| shows the dispersion curve along the path A 
(<7 y -direction) and B ((^-direction) near Q in Figs. || (a) 
and (b). The spin wave velocity v'£ pin (v^ pin ) is defined 
by the slope of the q x - (q y -) direction at u> ~ 0. The 
spin modulation parallel to the stripe (domain wall) cor- 
responds to v* pin . In this direction, the staggered spin 
moment has a constant amplitude. The modulation per- 
pendicular to the stripe corresponds to in . In this 
direction, the spin moment is modulated and suppressed 
when it crosses the stripe region. In Fig. ||, in > v^ pin . 
It indicates that the spin modulation is easier for the di- 
rection perpendicular to the stripe. In other words, the 
effective exchange integral J across the stripe becomes 
weaker than that parallel to the stripe. It is the first time 
to microscopically derive the anisotropy of the spin wave 
velocity. As U increases, u sp i n decreases. These results 
are reasonable in view of the correspondence between 
Hubbard and Heisenberg models: v sp i n ~ J oc t 2 /U. We 
have done the same calculation for the diagonal stripe to 
confirm that the spin wave velocity is similar. 

The longitudinal spin mode x'zz (<3j <£>) is shown in 
Fig. H along paths A and B of Fig. [l|. Low energy 
modes appears at (2n + 1)Q. Along path A, the dis- 
persion relation is continuous and repeated, touching 
at (2n + 1)Q. Away from Q, the intensity decreases 
as x£,(JQ,0.01t)/j&(Q,0.01t) - 5.9 x KT 2 (I = 3), 
6.9 x 1(T 4 (1 = 5), 5.3 x 1(T 4 (1 = 7). These ratios are 
larger than those given by Eq. (^), which is, thus, not sat- 
isfied for Xzz- The intensity x"*(q,w) ~ |x!L(q,w) near 
Q along each dispersion relation in our results. There, 
x" z (q, ~ 1/oj with increasing u>. 

The charge mode x'nn (Qj w ) i s shown in Fig. [| along 
paths C and D of Fig. fj] The low energy mode appears 
at 2nQ. It has the identical dispersion curve as that 
of Xzz, because S z and n couple each other in the RPA 
equation. As in Fig. ^ (a), it is a continuous curve along 
path C and its intensity deceases away from 2Q. The 2Q 
mode comes from even harmonics, as (n2Q 1 a) describes 
the charge modulation in the ISDW. 

We analyze the oscillation pattern of S z and n by Eq. 
(||). Along the dispersion curve of Xzz of Fig. || (xnn of 
Fig. U), the response is enhanced by the resonance with 
the external field coupled to S z (n). The response of 
5(S z (r,t)) and S(n(r,t)) shows large amplitude near the 
stripe region. It means that this excitation is related to 
the motion of the stripe, i.e., phason mode. This analysis 
shows that the collective mode at u> ~ is the transla- 
tional motion, where whole stripes move together to the 
same direction. When the pinning (such as the energy 
difference between the site-centered stripe and the bond- 
centered stripe) is negligible, this translational mode is 
a Goldstone mode and gapless as in Figs |] and [|. From 
the analysis of Eq. (^|), we understand that the excita- 
tion along the g^-direction is a compress mode and that 
along the ^-direction is a meandering mode. In the 
compress mode, the inter-stripe distance is modulated 



periodically in the direction perpendicular to the stripe, 
with keeping the straight line shape. In the meandering 
mode, each stripe meanders along the stripe direction 
with keeping the same inter-stripe distance. From Fig. 
||, we see i>p hason < "p hason - It means the meandering mo- 
tion is easier to occur compared with the compress mode 
in the vertical stripe case of this model. As N increases, 
^phason an ^ "phason decrease. It is because effective inter- 
action between neighbor stripes is weak and each stripe 
can move more freely when the inter-stripe distance be- 
comes long. This slow velocity of the longitudinal mode 
for large N may be related to the un-identified Fincher- 
Burke mode |13) observed in Cr. This identification de- 
serves further experimental and theoretical studies. 

We are also interested in the silent position q = S = 
2n(i ± 5, h), which is an equivalent position to Q in the 
paramagnetic state above Neel temperature (see Fig{l]). 
This silent mode is related to a critical scattering in 
Cr pl[] . There exists a large intensity in Xxx(&> whose 
dispersion has a gap of the order E g . In x" 2 (q, u;), the 
excitation at S (but slightly shifted to lower q x ) is shifted 
to lower energy asiV = 8^12^16. It almost touches 
u> = for N = 16. It may suggest that we are approach- 
ing the transition to another low energy ground state 
(such as diagonal stripe) for lower rih EllQ. 

So far we mention only the insulating stripe state. The 
corresponding metallic stripe state can be also stabilized 
by merely introducing the next nearest hopping t'. The 
lowest energy state is given by 5 ~ n^, and the Fermi 
level situates in the so-called midgap band @. Then, 
we set N = 8 in the metallic case. Since v^ pin increases 
and < P in decreases, we find v^ pin > v* pin in x" x . There is 
low energy excitation also at S. In %" z and x'nm the exci- 
tation at Q or 2Q has a gap. The low energy excitation 
appears along the line at q x = tt/2 which is presented 
by a line in Fig. |. It is the ID CDW or SDW fluctua- 
tion mode within the stripe line, and originated from the 
Fermi surface nesting 2kpiD of the parallel ID Fermi lines 
(Fig. 12(d) in Ref. @) of the stripe state |]|. Since 
the ID Fermi state has a gap near (f, f ), the intensity of 
x"(2kpiD, w ~ 0) vanishes near (f , 7r). These low energy 
excitations are diffusive since E g — in metallic state. 

In summary, we have investigated the dynamical sus- 
ceptibilities of transverse and longitudinal spin channels 
and charge one for the whole space spanned by 2D wave 
vector q and the energy lo, and identified several ele- 
mentary excitations; the spin wave mode and the phason 
mode related to the motion of the stripe line. This al- 
lows us to construct the whole landscape of (q, u>) space 
for the excitation spectra of various channels: The iden- 
tical dispersion relation is replicated at every (2n+ 1)Q, 
which has the anisotropic excitation cones along q x - and 
Qy-directions. Our predictions about x(q, u>) are directly 
testable by careful inelastic neutron experiments on Cr 
alloys and underdoped cuprates. 
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FIG. 2. Spin transverse mode Xx X {<i,oj) along path A (a) 
and B (b). u) is scaled by t. We cut off the peak height for 
Xxx > 10 to show the low intensity structure. 




FIG. 3. Dispersion curve of the spin-wave excitation (vg pin 
for path B and v^ pin for path A) and phason excitation (Vp hason 
for path B and «p hason for path A). The slope at lu ~ gives 
velocity of each mode. 
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FIG. 1. The paths A-D along which we show in 
the momentum space. The point Q is the ordering vector. 
2Q, 3Q, ■ ■ ■ are its higher harmonics points. S is the silent 
position. The line 21cfid shows the nesting wave number of 
the ID Fermi surface in the metallic stripe state. 



FIG. 4. Spin longitudinal mode x"z( c l>' i; ) along path A 
(a) and B (b). We cut off the peak height for x'L > 200. 
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FIG. 5. Charge mode Xnni^l^) along path C (a) and D 



(b). 
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